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We present a measurement of the top-quark width in the lepton + jets decay channel of ti events 
produced in pp collisions at Fermilab's Tevatron collider and collected by the CDF II detector. From 
a data sample corresponding to 4.3 fb -1 of integrated luminosity, we identify 756 candidate events. 
The top-quark mass and the mass of the hadronically decaying W boson that comes from the top- 
quark decay are reconstructed for each event and compared with templates of different top-quark 
widths (Ft) and deviations from nominal jet energy scale (Ajes) to perform a simultaneous fit for 
both parameters, where Ajes is used for the in situ calibration of the jet energy scale. By applying 
a Feldman-Cousins approach, we establish an upper limit at 95% confidence level (CL) of V t < 7.6 
GeV and a two-sided 68% CL interval of 0.3 GeV < F t < 4.4 GeV for a top-quark mass of 172.5 
GeV/c 2 , which are consistent with the standard model prediction. 

PACS numbers: 14.65.Ha, 13.85.Qk, 12.15.Ff 
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The top quark is the heaviest known elementary par- 
ticle, whose large mass results in the largest decay width 
and hence the shortest lifetime of the quarks in the stan- 
dard model (SM). A precise measurement of the top- 
quark width r t is a good test of the standard model, 
whose prediction at the Born level [l| is affected by 
the Quantum Chromodynamics (QCD) radiative correc- 
tions of order 10% [2j, as well as by electroweak correc- 
tions [HIH, which are of order 1.5%. The dominant decay 
mode of the top quark in the SM produces a W boson 
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and a bottom quark (6). At leading order the total top- 
quark width is given by T° = \V tb \ 2 G F m 3 t /(8irV2), where 
Vtb, Gf and m t are the Cabibbo-Kobayashi-Maskawa 
matrix element, Fermi coupling constant and the top- 
quark mass, respectively. If we takelV^I to be unity, the 
next-to-leading order calculation [l], [2| with QCD and 
electroweak corrections predicts T t of 1.3 GeV at a top- 
quark mass of 172.5 GeV/c 2 Q and this corresponds to 
a lifetime of 5 x 1CP 25 s. 

A deviation from the SM could indicate a significant 
contribution of non-SM particles. Novel top-quark de- 
cay modes motivated by the large top-quark mass in- 
clude decay to a charged Higgs t — > b + H + , decay 
to its supersymmetric scalar partner stop plus neutrali- 
nos [13, El) an d flavor changing neutral current (FCNC) 
top-quark decays [HI]. Therefore, the direct measure- 
ment of r t is a general way to constrain such processes. 
The first direct measurement of T t was carried out with 
an integrated luminosit y o f 1 fb _1 of CDF data in the 
lcpton + jets channel |l3| and set an upper limit on 
r t < 13.1 GeV at 95% confidence level (CL), while the 
result of a recent analysis from the DO experiment at the 
Tevatron quotes an indirect top-quark width measure- 
ment of T t = 1.99±g:|| GeV [I| . In this report of the 
second direct measurement of T t , we increase the CDF 
data set to 4.3 fb _1 in the lepton + jets channel, apply 
a kernel density estimation (KDE) technique [HI, [Tfl] to 
make templates, determine the jet energy scale (JES) cal- 
ibration in situ, and use new methods for setting and in- 
corporating systematic effects. We set a two-sided bound 
on the top-quark width at 68% CL for the first time. 

CDF II [ItJ is a general-purpose detector located at 
one of the two collision points along the ring of the 
Tevatron accelerator. A silicon microstrip tracker and a 
cylindrical drift chamber in a 1.4 T magnetic field serve 
as a charged particle tracking system. Electromagnetic 
and hadronic calorimeters are used to measure the en- 
ergies of electrons and jets. Outside the calorimeters lie 
drift chambers which can detect muons. We employ a 
cylindrical coordinate system for the detector where 8 
and <j) are the polar and azimuthal angles, respectively, 
with respect to the proton beam, and pseudorapidity 
rj = — lntan(0/2). Transverse energy and momentum 
are defined as Et — Esm9 and pt = psin#, respec- 
tively, where E and p are energy and momentum. 

Top quarks decay almost exclusively to a W boson and 
a b quark through the weak interaction in the SM. We 
identify tt events in the lepton + jets channel, where 
one W boson decays to a charged lepton and neutrino, 
and the other W boson decays to two quarks. The tt 
candidate events used in this analysis are collected by 
triggers that identify at least one high-p^ lepton. Of- 
fline these events are selected by requiring a high-Er 
electron or high-pT muon (Et or px > 20 GeV), large 
missing transverse energy #t [£t > 20 GeV) due to 
the undetected neutrino from the leptonic W decay, and 
at least four hadronic jets. Jets are reconstructed with 
the jetclu [3| cone algorithm using a cone radius of 



AR = y/Aif + Ac/} 2 = 0.4. To determine if a jet comes 
from a b quark, the SECVTX [19| algorithm, which makes 
use of the transverse decay length of a & quark inside a 
jet (6-tag), is applied. At least one jet must be iden- 
tified as &-tagged. We divide the candidate events into 
those with one &-tagged jet and those with two or more 
6-tagged jets in order to improve the usage of statistical 
information, since these two kinds of events have differ- 
ent signal-to-background ratios. When an event has one 
&-tagged jet (6-jet), we require this event to have exactly 
four jets each with Et > 20 GeV; when an event contains 
two or more &-jets, three jets are required to have Et > 
20 GeV, the fourth must have E T > 12 GeV, and the 
event is allowed to have extra jets. More details about 
event selection criteria can be found in (20| . 

Monte Carlo (MC) simulated signal samples are cre- 
ated for a fixed top-quark mass of 172.5 GeV/c 2 by the 
PYTHIA version 6.216 [2l| event generator and have dif- 
ferent values of T t between 0.1 GeV and 30 GeV, as well 
as various values of Aj^s, which is the difference be- 
tween the JES effects in MC simulation and data and 
has a range from -3.0 a c to +3.0 a c , where a c is the 
CDF JES fractional uncertainty [13] • The overall rate of 
background events with one W boson and additional jets 
(W + jets), the dominant background process, is deter- 
mined using data after subtracting off the rate of events 
coming from QCD multi-jet production (non-W events), 
and separating out a MC based estimate for electroweak 
processes (EWK) such as diboson (WW/WZ/ZZ) and 
single-top production. The fractions of W + jets events 
with heavy flavor quarks (Wc, Wcc, and Wbb events) are 
determined from MC simulated samples. The rate with 
which events with a W boson and light flavor quarks con- 
tain a misidentified 6-jet is determined using data sam- 
ples triggered by the presence of jets. Table U summa- 
rizes the background compositions, and the selection cri- 
teria for determining the background rates are described 
in |23j. Diboson backgrounds are modeled with PYTHIA 
version 6.216 [2l| and W + jets by ALPGEN version 
2.10 (24|, with jet fragmentation modeled by PYTHIA ver- 
sion 6.325 [2l|. Sin gle- top production events are gener- 
ated by MADEVENT [25|] and their fragmentation is mod- 
eled with pythia version 6.409 [26j . 

We use a template method to extract T t . Two ob- 
servables, the reconstructed top-quark mass (m[ eco ) and 
the invariant mass of the two jets from the hadronically 
decaying W boson (rrijj), are built for each data event 
or MC simulated event (both signal and background). 
With the assumption that the leading (highest Et) four 
jets in the detector come from the four primary quarks of 
tt events in lepton + jets channel, there are 12 possible 
assignments of jets to quarks in each event. The neutrino 
transverse momentum is calculated from the imbalance 
of the transverse momentum of decaying products, jets 
and lepton, with unclustered energy taken into account, 
which is the energy in the calorimeter not associated with 
the lepton or one of the four leading jets. We use a x 2 -like 
kinematic fitter [27] to fit the top-quark mass for each as- 



TABLE I: The sources and expected numbers of background 
events in the lepton + jets channel, and the number of events 
observed for single 6-tag and double b-tag samples after event 
selection, x 2 cut, and boundary cuts. 



TABLE II: Summary of changes in measured Ft due to sys- 
tematic effects. 





single &-tag 


double 6-tag 


W + jets 


85.6 ± 21.8 


9.8 ± 2.9 


non-IV 


24.5 ± 20.6 


2.4 ± 1.8 


EWK 


10.2 ± 0.8 


2.4 ± 0.2 


Total background 
Observed events 


120.2 ± 30.0 
542 


14.6 ± 3.4 
214 





— I> 1.5 GeV 




, — I>10.0GeV 




] I>20.0 GeV 




1 (a) 







JO 0.02 

o 

CL. 



150 200 250 300 350 

m{ eco (GeV/c 2 ) 





--A JES = -2.0 o c 








\ \ . %s=2-»». 




\V\ (b) 



m (GeV/c ) 



FIG. 1: (a) Probability density functions of m T t ec ° from dou- 
ble b-tag events for MC simulated samples of different values 
of Ft; (b) p.d.f.'s of rrijj from double b-tag events for MC 
simulated samples of different values of Ajes. 



signment, assuming the mass equality of the top and anti- 
top quarks, and take m r t eco from the assignment that has 
the lowest \ 2 - Events with x 2 > 9.0 are removed from the 
sample to reject poorly reconstructed events. We also ap- 
ply boundary cuts on m r t eco (110 GeV/c 2 < m r t eco < 350 
GeV/c 2 ) and rrijj (50 GeV/c 2 < m n < 115 GeV/c 2 for 
single 6-tag events and 50 GeV/c 2 < rrijj < 125 GeV/c 2 
for double 6-tag events) and normalize the probability 
density functions (p.d.f.) in these regions. The di-jet 
mass rrijj is calculated as the invariant mass of two non- 
6-tagged jets which provides the closest value to the world 
average W boson mass of 80.40 GeV/c 2 [28| . The esti- 
mated number of background events and observed num- 
ber of events from a data set corresponding to an in- 
tegrated luminosity of 4.3 fb _1 after event selection, % 2 
cut, and boundary cuts are listed in Table fl] After event 
reconstruction, we use the MC simulated models of sig- 
nal and background processes to build two-dimensional 
p.d.f.'s that give the probabilities of observing a pair of 
values of m r t eco and rrijj, given T t and Ajes- We em- 
ploy a KDE that associates to each data point a func- 
tion (called a kernel function) and uses a non-parametric 
method to estimate the p.d.f.'s of a variable by summing 
all the kernel functions, without any assumption about 
the functional form of the p.d.f.'s. Figure Q] shows the 
p.d.f.'s of m T t eco with different T t and the rrijj with var- 
ious Ajes from a full MC simulation. We compare the 



Systematic Sources 


Ar top (GeV) 


Jet energy resolution 


1.1 


Color Reconnection 


0.9 


Generator 


0.4 


Residual JES 


0.3 


Parton distribution functions 


0.3 


Multiple Hadron Interaction 


0.3 


gluon gluon fraction 


0.3 


Initial/Final state radiation 


0.2 


Lepton energy scale 


0.2 


6 jet energy 


0.2 


Background shape 


0.1 


Total systematic effect 


1.6 



distributions of data with signal and background p.d.f.'s 
using an unbinned maximum likelihood fit [29| . where 
the likelihood function C is the same as what is defined 
in [lj|] except that the fitting parameters are now T t and 
Ajes- We minimize the negative logarithm of the like- 
lihood using MINUIT [3(| to extract the top-quark width. 
The fitting to Ajes reduces the JES systematic effect on 
Tt and thus improves the sensitivity to the top width. 

We set the limit(s) on T t via the Feldman- Cousins 
method 31( which determines the confidence intervals. 
The ordering parameter for MC simulated samples that 
appears in [3jj| is defined here as Ax 2 = xlnput ~ Xmin> 
where x 2 = —2 log(£) (different from the x 2 mentioned in 
event reconstruction), Xmin ^ s ^ ne minimal x 2 value and 
Xinput 1S the X 2 a t the real value of parameters T t and 
Ajes of the MC simulated sample. We project the likeli- 
hood function C onto the T t axis [32| . For each value of T t 
we run 6000 pseudo-experiments that generate a distribu- 
tion of Ax 2 from which we calculate a critical value Ax 2 
so that 95% of the pseudo-experiments have a A% 2 falling 
in the interval [0, A% 2 ]- With MC simulated samples of 
21 different top widths T t we get a profile of Ax 2 (r t ). 
When analyzing the data we obtain Ax 2 (r t |data) = 
— 21og(£) + 21og(£ ), where Cq is the maximum like- 
lihood value of data fitting, then Ax 2 (r t |data) is com- 
pared with Ax 2 (rt) and the accepted interval of Tt is 
all points such that Ax 2 (r t |data) < Ax 2 (r t ). From the 
above method we obtain a purely statistical upper limit 
on T t at 95% CL, T t < 6.7 GeV and a two-sided limit of 
0.5 GeV < T t < 3.9 GeV at 68% CL. 

We examine systematic effects by comparing MC sim- 
ulated experiments in which we float parameters within 
their uncertainties. As seen from Table ILTl the dominant 
systematic effects come from jet energy resolution and 
color reconnection (CR) [H, [34[ , which is a rearrange- 
ment of the underlying color structure of an event from 
its simplest configuration. For the jet energy resolution 
effect, we compare jet energy resolution between data 
and MC simulated samples using one photon + one jet 
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CDF II (4.3 fb" 1 ) 
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FIG. 2: (a) Contours of the two-dimensional negative log 
likelihood function from data fit. The three different con- 
tours represent different values of — log(£): 0.5, 2.0, and 4.5. 
Systematic effects are not included here, (b) Overlap of the 
Axc(r«) profile and the data fit that comes from projection 
of the two-dimensional data fit onto the Tt axis, the intersec- 
tion of which gives a limit (s) on P t . Systematic effects are 
included in the plots, both for 68% and 95% CL. 

events and smear jet energy with the difference between 
data and MC simulated samples. We study the effect 
of CR by using pythia version 6.4 with different tunes 
(with and without CR) and evaluate the difference. The 
systematic effect due to JES is very small because we 
perform an in situ JES calibration. Other smaller sys- 
tematic effects include those due to the MC generator, 
the parton distribution functions, and multiple hadron 
interactions, details of which can be found in [j| l35j . 
The total change of measured T t due to these systematic 
effects is 1.6 GeV. We studied the dominant systematic 
uncertainties by varying top-quark width, and found no 
significant dependence of systematic effects on different 



top-quark widths. 

To incorporate systematic effects into the limit (s) on 
r t we use a convolution method for folding systematic 
effects into the likelihood function [3(1 H?} ■ We convolve 
the likelihood function with a Gaussian p.d.f. that has a 
width equal to 1.6 GeV and is centered at 0. With this 
new likelihood function we apply the Feldman- Cousins 
approach and find an upper limit of T t < 7.6 GeV at 
95% CL. Using the same approach we are also able to 
set a two-sided bound for T t at 68% CL: 0.3 GeV < 
T t < 4.4 GeV. Figure HJa) shows the data fit from the 
two-dimensional likelihood function with the statistical 
uncertainty. The overlap of the Ax 2 (r t ) profile and the 
one-dimensional data fit that comes from the projection 
of the two-dimensional likelihood function is shown in 
Fig. Hfb), on which the point (s) of interception gives the 
limit (s) of T t . 

In conclusion, a top-quark width measurement in the 
lepton + jets channel is presented. Using a data set cor- 
responding to an integrated luminosity of 4.3 fb -1 col- 
lected by CDF and an in situ JES calibration, we set an 
upper limit T t < 7.6 GeV at 95% CL assuming a top- 
quark mass M top — 172.5 GeV/c 2 , which is consistent 
with the standard model. We also quote 0.3 GeV < r t < 

4.4 GeV at 68% CL, which corresponds to a life time of 

1.5 xl0~ 25 s < n < 2.2 xlO -24 s. For a typical quark 
hadronization time scale of 3.3 xl0~ 24 s (corresponding 
to 200 MeV) [3Sll39j . our result supports top-quark decay 
before hadronization. 
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